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Summary. — An overview of the ab initio no-core shell model is presented. Recent 
results for light nuclei obtained with the chiral two-nucleon and three-nucleon inter- 
actions are highlighted. Cross section calculations of capture reactions important 
for astrophysics are discussed. The extension of the ab initio no-core shell model 
to the description of nuclear reactions by the resonating group method technique is 
outlined. 



1. Introduction 

The major outstanding problem in nuclear physics is to calculate properties of finite 
nuclei starting from the basic interactions among nucleons. This problem has two parts. 
First, the basic interactions among nucleons are complicated, they are not uniquely de- 
fined and there is evidence that more than just two-nucleon forces are important. Second, 
the nuclear many-body problem is very difficult to solve. This is a direct consequence 
of the complex nature of the inter- nucleon interactions. Both short-range and medium- 
range correlations among nucleons are important and for some observables long-range 
correlations also play a significant role. Various methods have been used to solve the 
few- nucleon problem in the past. The Faddeev method IJ has been successfully applied 
to solve the three-nucleon bound-state problem for different nucleon-nucleon potentials 
[2l[3l|4] . For the solution of the four-nucleon problem one can employ Yakubovsky's gener- 
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alization of the Faddeev formalism [5] as done, e.g., in Refs. [6^ or f7j. Alternatively, other 
methods have also been succesfully used, such as, the correlated hypcrspherical harmon- 
ics expansion method [8l[9] or the Green's function Monte Carlo method [10]. Recently, 
a benchmark calculation by seven different methods was performed for a four-nucleon 
bound state problem [11]. However, there are few approaches that can be successfully 
applied to systems of more than four nucleons when realistic inter-nucleon interactions 
are used. Apart from the coupled cluster method [HI [T31 [T?] [IS] applicable typically to 
closed-shell nuclei, the Green's function Monte Carlo method is a prominent approach 
capable to solve the nuclear many-body problem with realistic interactions for systems 
of up to A = 12. Another method developed recently applicable to light nuclei up to 
A = 16 and beyond is the ab initio no-core shell model (NCSM) [TOj. In this paper, an 
overview of this approach is given and results obtained very recently are presented. Also, 
future developments, in particular applications to nuclear reactions, are outlined. 



2. ~ Ab initio no-core shell model 

In the ah initio no-core shell model, we consider a system of A point-like non- 
relativistic nucleons that interact by realistic two- or two- plus three-nucleon (NNN) 
interactions. Under the term realistic two-nucleon (NN) interactions we mean NN po- 
tentials that fit nucleon-nucleon phase shifts with high precision up to certain energy, 
typically up to 350 MeV. A realistic NNN interaction includes terms related to two-pion 
exchanges with an intermediate delta excitation. In the NCSM, all the nucleons are 
considered active, there is no inert core like in standard shell model calculations. There- 
fore the "no-core" in the name of the approach. There are two other major features in 
addition to the employment of realistic NN or NN-I-NNN interactions. The first one is 
the use of the harmonic oscillator (HO) basis truncated by a chosen maximal total HO 
energy of the A-nucleon system. The reason behind the choice of the HO basis is the fact 
that this is the only basis that allows to use single-nucleon coordinates and consequently 
the second-quantization representation without violating the translational invariance of 
the system. The powerful techniques based on the second quantization and developed for 
standard shell model calculations can then be utilized. Therefore the "shell model" in the 
name of the approach. As a downside, one has to face the consequences of the incorrect 
asymptotic behavior of the HO basis. The second feature comes as a consequence of the 
basis truncation. In order to speed up convergence with the basis enlargement, we con- 
struct an effective interaction from the original realistic NN or NN-I-NNN potentials by 
means of a unitary transformation. The effective interaction depends on the basis trun- 
cation and by construction becomes the original realistic NN or NN-I-NNN interaction as 
the size of the basis approaches infinity. In principle, one can also perform calculations 
with the unmodified, "bare", original interactions. Such calculations are then variational 
with the HO frequency and the basis truncation parameter as variational parameters. 
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The starting Hamiltonian of the ab initio NCSM is 
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where m is the nucleon mass, l^N.y the NN interaction, VNNN.ijfc the three-nucleon 
interaction. In the NCSM, we employ a large but finite HO basis. Due to properties of the 
realistic nuclear interaction in Eq. ([1]), we must derive an effective interaction appropriate 
for the basis truncation. To facilitate the derivation of the effective interaction, we 
modify the Hamiltonian ([1]) by adding to it the centcr-of-mass (CM) HO Hamiltonian 
HcM = TcM + C/cM, where Ucm = ^Amn^lP, R=\ Y.t=i The effect of the HO CM 
Hamiltonian will later be subtracted out in the final many-body calculation. Due to the 
translational invariance of the Hamiltonian ([1]) the HO CM Hamiltonian has in fact no 
effect on the intrinsic properties of the system. The modified Hamiltonian can be cast 
into the form 



A r 
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2'2. Basis. - In the ab initio NCSM, we use a HO basis. A single-nucleon HO wave 
function can be written as 



(3) 



'Pnlm{r;b) = Rnl{r\b)Yira{f) 



with Rni {r, b) the radial HO wave function and b the HO length parameter related to the 

HO frequency as = with m the nucleon mass. The HO length parameter b is 

often dropped in Rni and (pnim in the following text to simplify notation. 

The HO wave functions have important transformation properties that we utilize 
frequently. Let us consider two particles with different masses of ratio d = m2/mi 
moving in a HO well. Their relative and center-of-mass coordinates can be defined by 
an orthogonal transformation 



(4) 
(5) 



R = 
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7^1 



1 



-J2 , 



1 + d 
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where all the vectors in ^ and (5) are defined as products of square root of the respective 
mass and the position vector: f = ^/mx. The product of the single-particle HO wave 
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functions can then be expressed as a linear combination of the relative-coordinate HO 
wave function and the CM coordinate HO wave function: 



(6) ['/'«i/i(ri)</'«2/2(^2)]i^' = 



E 

alNL 



{nlNLK\n^hn2l2K)d['Pni{^^NL{R)]f^ 



with {nlN LK\nilin2l2K) d a generalized HO bracket that can be evaluated according to 
the algorithm in, e.g. Ref. [T7] . 

As the NN and NNN interactions depend on relative coordinates and/or momenta, 
the natural coordinates in the nuclear problem are the relative, or Jacobi, coordinates. 

We work in the isospin formalism and consider nucleons with the mass m. A gener- 
alization to the proton-neutron formalism with unequal masses for the proton and the 
neutron is straightforward. We will use Jacobi coordinates that are introduced as an or- 
thogonal transformation of the single-nucleon coordinates. In general, Jacobi coordinates 
are proportional to differences of centers of mass of nucleon sub-clusters. 

For the present purposes we consider just a single set of Jacobi coordinates. More 
general discussion can be found in Ref. [THj. The following set 



(7) 
(8) 
(9) 

(10) 
(11) 



[ri + r2 + . . . + rA] , 




is useful for the construction of the antisymmetrized HO basis. Here, is proportional 
to the center of mass of the A-nucleon system. On the other hand, is proportional to 
the relative position of the p + 1-st nucleon and the center of mass of the p nucleons. 

As nucleons are fermions, we need to construct an antisymmetrized basis. Here we 
illustrate how to do this for the simplest case of three nucleons. One starts by introducing 
a HO basis that depends on Jacobi coordinates and ^2, defined in Eqs. (8) and (9), 
e.g., 



(12) 



\{nlsjt-MCJ)JT) 



Here n, I and A/", C are the HO quantum numbers corresponding to the harmonic os- 
cillators associated with the coordinates (and the corresponding momenta) and ^2, 
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respectively. The quantum numbers s,t,j describe the spin, isospin and angular mo- 
mentum of the relative-coordinate two-nucleon channel of nucleons 1 and 2, while is 
the angular momentum of the third nucleon relative to the center of mass of nucleons 1 
and 2. The J and T are the total angular momentum and the total isospin, respectively. 
Note that the basis (fT2|) is antisymmetrized with respect to the exchanges of nucleons 
1 and 2, as the two-nucleon channel quantum numbers are restricted by the condition 
(— 1)'+'*+* = —1. It is not, however, antisymmetrized with respect to the exchanges of 
nucleons 1 <-> 3 and 2 <-> 3. In order to construct a completely antisymmetrized basis, 
one needs to obtain eigenvectors of the antisymmetrizer 

(13) A' = i(l + T(-)+TW), 

where T(+' and T*^~' are the cyclic and the anti-cyclic permutation operators, respec- 
tively. The antisymmetrizer X is a. projector satisfying XX = X. When diagonalized 
in the basis (fT2|) . its eigenvectors span two eigenspaces. One, corresponding to the 
eigenvalue 1, is formed by physical, completely antisymmetrized states and the other, 
corresponding to the eigenvalue 0, is formed by spurious states. There are about twice 
as many spurious states as the physical ones [19]. 

Due to the antisymmetry with respect to the exchanges 1 <-> 2, the matrix elements in 
the basis (fT^ of the antisymmetrizer X can be evaluated simply as {X) = ^(1 — 2(P2,3))7 
where P2,3 is the transposition operator corresponding to the exchange of nucleons 2 and 
3. Its matrix element can be evaluated in a straightforward way, e.g.. 



{inillSijih-N'iCiJi)JT\P2,3\{n2l2S2j2t2;J^2C2j2)JT) 

1 1 



2 



T t2 
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h Si ji ] C ^2 S2 j2 

xY,L^S^JlhJlJ2hS2i-l)'' \ Ci i Jl \ \ C2 i J2 

LS [lsj][lsj 

1 1 



(14) x\ I I j> {nihJ\fiCiL\U2C2n2l2L)3 



where Ni = 2ni + h + 2M+Ci,i = 1,2; j = y/2j + 1; and {nihNiCiL\M2C2n2l2L)z is the 
general HO bracket for two particles with mass ratio 3 as defined, e.g., in Ref. [17]. The 
expression pi)) can be derived by examining the action of P2,3 on the basis states ([T^ . 
That operator changes the state |nZ(fi), A/'£(f2), L) to \nl{£^\),N C{^' L) ^ where f'^, i — 
1,2 are defined as = 1,2 but with the single-nucleon indexes 2 and 3 exchanged. 
The primed Jacobi coordinates can be expressed as an orthogonal transformation of 
the unprimed ones, see Eq. ([6]). Consequently, the HO wave functions depending on 
the primed Jacobi coordinates can be expressed as an orthogonal transformation of the 
original HO wave functions. Elements of the transformation are the generalized HO 
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brackets for two particles with the mass ratio d, with d determined from the orthogonal 
transformation of the coordinates, see Eq. (jl]). 

The resulting antisymmetrized states can be classified and expanded in terms of the 
original basis as follows 

(15) \NiJT) = ^{nlsjt;JVCJ\\NiJT)\{nlsjt;J\fCJ)JT) , 

where N = 2n+l+2Af+C and where we introduced an additional quantum number i that 
distinguishes states with the same set of quantum numbers N, J, T, e.g., i ~ 1, 2, ... r 
with r the total number of antisymmetrized states for a given J, T . The symbol 
{nlsjt\N CJ\\NiJT) is a coefficient of fractional parentage. 

A generalization to systems of more than three nucleons can be done as shown e.g. 
in Ref. [18] . It is obvious, however, that as we increase the number of nucleons, the 
antisymmetrization becomes more and more involved. Consequently, in the standard 
shell model calculations one utilizes antisymmetrized wave functions constructed in a 
straightforward way as Slater determinants of single-nucleon wave functions depending 
on single-nucleon coordinates ipi{fi). It follows from the transformations ([6]) that a use 
of a Slater determinant basis constructed from single nucleon HO wave functions such as 

(16) 

results in eigenstates of a translationally invariant Hamiltonian that factorize as products 
of a wave function depending on relative coordinates and a wave function depending on 
the CM coordinates. This is true as long as the basis truncation is done by a chosen 
maximum of the sum of all HO excitations, i.e.: X^iLil-^^-i + h) < Ntotmax- In Eq- (HH), 
a and t are spin and isospin coordinates of the nucleon. The physical eigenstates of an 
translationally invariant Hamiltonian can then be selected as eigenstates with the CM in 
the Oh^l state: 

(fi . . . fAcri ■ . . a- ATI . . . ta\AXJMTMt)su 
(17) = (fi . . . Ca-iCTi . . . <JATi . . . ta\AXJMTMt)^ooo{^o; h) , 

For a general single-nucleon wave function this factorization is not possible. A use of any 
other single nucleon wave function than the HO wave function will result in mixing of 
CM and internal motion. 

In the ah initio NCSM calculations, we use both the Jacobi-coordinate HO basis and 
the single-nucleon Slater determinant HO basis. One can choose, whichever is more 
convenient for the problem to be solved. One can also mix the two types of bases. In 
general for systems of A < 4, the Jacobi coordinate basis is more efficient as one can 
perform the antisymmetrization easily. The CM degrees of freedom can be explicitly 
removed and a coupled J^T basis can be utilized with matrix dimensions of the order 
of thousands. For systems with A > 4, it is in general more efficient to use the Slater 
determinant HO basis. In fact, we use so-called m-scheme basis with conserved quantum 
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numbers M — X]i=i "^ii parity tt and Mt = X]i=i "^tj- The antisymmetrization is trivial, 
but the dimensions can be huge as the CM degrees of freedom are present and no JT 
coupUng is considered. The advantage is the possibihty to utihze the powerful second 
quantization technique, shell model codes, transition density codes and so on. 

As mentioned above, the model space truncation is always done using the condition 
Z]iLi(2"-j + li) < Ntotmax- Often, instead of Ntotmax we introduce the parameter A'^max 
that measures the maximal allowed HO excitation energy above the unperturbed ground 
state. For A = 3,4 systems iVmax — Ntotmax- For the p-shell nuclei they differ, e.g. for 

^Li, Amax — Ntotmax — 2, for "'^^C, A^max — Ntotmax ~ 8 CtC. 

2'3. Effective interaction. - In the ah initio NCSM calculations we use a truncated HO 
basis as discussed in previous sections. The inter-nucleon interactions act, however, in the 
full space. In order to obtain meaningful results in the truncated space, or model space, 
the inter-nucleon interactions needs to be renormalized. We need to construct an effective 
Hamiltonian with the inter-nucleon interactions replaced by effective interactions. By 
meaningful results we understand results as close as possible to the full space exact results 
for a subset of eigenstates. Mathematically we can construct an effective Hamiltonian 
that exactly reproduces the full space results for a subset of eigenstates. In practice, we 
cannot in general construct this exact effective Hamiltonian for the A-nucleon problem 
we want to solve. However, we can construct an effective Hamiltonian that is exact for 
a two-nucleon system or for a three-nucleon system or even for a four-nucleon system. 
The corresponding effective interactions can then be used in the A-nucleon calculations. 
Their use in general improves the convergence of the problem to the exact full space 
result with the increase of the basis size. By construction, these effective interactions 
converge to the full-space inter-nucleon interactions therefore guaranteeing convergence 
to exact solution when the basis size approaches the infinite full space. 

In our approach we employ the so-called Lee-Suzuki similarity transformation method 
PPI I21j . which yields a starting-energy independent hermitian effective interaction. We 
first recapitulate general formulation and basic results of this method. Applications of 
this method for computation of two- or three-body effective interactions are described 
afterwards. 

Let us consider an arbitrary Hamiltonian H with the eigensystem Ek, \k), i.e., 



Let us further divide the full space into the model space defined by a projector P and the 
complementary space defined by a projector Q, P+Q = 1. A similarity transformation of 
the Hamiltonian He'^ can be introduced with a transformation operator w satisfying 
the condition uj = QutP. The transformation operator is then determined from the 
requirement of decoupling of the Q-space and the model space as follows 



(18) 



H\k)^Ek\k) . 



(19) 



Qe^'^He'^P . 
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If we denote the model space basis states as \ap), and those which belong to the Q-space, 
as lag), then the relation Qe^'^ He'^ P\k) — 0, following from Eq. (fT9|) . will be satisfied 
for a particular eigenvector |fc) of the Hamiltonian (jlSp . if its Q-space components can be 
expressed as a combination of its P-space components with the help of the transformation 
operator uj, i.e., 

(20) (aglfc) =^(aQ|c^|ap)(ap|fc) . 

a p 

If the dimension of the model space is dp, we may choose a set K, of dp eigenevectors, 
for which the relation (|20p will be satisfied. Under the condition that the dp x dp 
matrix defined by matrix elements {ap\k) for |fc) G /C is invertible, the operator lu can 
be determined from (j20|) as 

(21) {aQ\uj\ap) = Y,{^Q\k)Ck\ap) , 

keK 

where we denote by tilde the inverted matrix of (ap|fc), e.g., J2ap(k\'^p)i'^p\k') ~ ^k,k' ^ 
for fc, k' g /C. 

The hermitian effective Hamiltonian defined on the model space P is then given by 

m 

(22) iJeff = [P(l + uj'^uj)P\ ^'"^ PHiP + QujP) [P(1 + Juj)P] '^'^ . 

By making use of the properties of the operator w, the effective Hamiltonian /?eff can be 
rewritten in an explicitly hermitian form as 

(23) Heff - [P{1 + uo^uj)P] '^'^ {P + Pu^Q)H{QujP + P) [P{1 + LJ^U0)P] '^'^ . 

With the help of the solution for lu (|2ip we obtain a simple expression for the matrix 
elements of the effective Hamiltonian 

k^K otpff a pttt 

(24) x{ap>»\{l+u)^ujy^''^\ap:) . 

For computation of the matrix elements of (1 + cj^w)^^/^, we can use the relation 

(25) {ap\{l+Lo^uj)\apn) = ^(ap|fc)(fc|ap,.) , 

to remove the summation over the Q-space basis states. The effective Hamiltonian ([M]) 
reproduces the eigenenergies Ek, k € IC in the model space. 



Ab initio NO-CORE SHELL MODEL CALCULATIONS FOR LIGHT NUCLEI 



9 



It has been shown [22] that the hermitian effective Hamiltonian ((23|) can be obtained 
directly by a unitary transformation of the original Hamiltonian: 

(26) i?off = Pe-'^He^P , 

with an anti- hermitian operator S = arctanh(aj — lo'^). The transformed Hamiltonian 
then satisfies decoupling conditions Qe^^ He^ P = Pe^^ He^ Q — 0. 

We can see from Eqs. (|24|) and (|25|) that in order to construct the effective Hamiltonian 
we need to know a subset of exact eigenvalues and model space projections of a subset 
of exact eigenvectors. This may suggest that the method is rather impractical. Also, it 
follows from Eq. ((24)) that the effective Hamiltonian contains many-body terms, in fact 
for an ^-nucleon system all terms up to A-body will in general appear in the effective 
Hamiltonian even if the original Hamiltonian consisted of just two-body or two- plus 
three-body terms. 

In the ab initio NCSM we use the above effective interaction theory as follows. Since 
the two-body part dominates the ^-nucleon Hamiltonian ([21), it is reasonable to expect 
that a two-body effective interaction that takes into account full space two-nuclcon cor- 
relations would be the most important part of the exact effective interaction. If the 
NNN interaction is taken into account, a three-body effective interaction that takes into 
account full space three-nucleon correlations would be a good approximation to the exact 
v4-body effective interaction. We construct the two-body or three-body effective inter- 
action by application of the above described Lee-Suzuki procedure to a two-nucleon or 
three-nucleon system. The resulting effective interaction is then exact for the two- or 
three-nucleon system. It is an approximation of the exact A-nucleon effective interaction. 

Using the notation of Eq.([2|), the two-nuclcon effective interaction is obtained as 

(27) y2cff,i2 = P2[e'^''{hi + h2 + l/i'^'^)e^^^ - (/ii + h2)]P2 , 

with 5*12 — arctanh(Li;i2 — i^\2) and P2 is a two-nucleon model space projector. The 
two-nucleon model space is defined by a truncation iVi2max corresponding to the A- 
nucleon iVmax- For example, for ^ = 3,4, iVi2max — Mnax, for p-shell nuclei with ^ > 5 
.^i2max = .^max + 2. Thc Operator 1^12 is obtained with the help of Eq. (|2ip from 
exact solutions of the Hamiltonian hi + h2 + which are straightforward to find. In 
practice, we actually do not need to calculate (jJi2, rather we apply Eqs. (|24p and ([^5]) 
with the two-nucleon solutions to directly calculate P2e^^^^{hi + h2 + V^'^)e^'-^P2- To 
be explicit, the two-nucleon calculation is done with 

(28) H'^ = H02 + V;^--" = ^ + Un'r' + VMN{V2r-) - , 

2m 2 A 

where r — \J\{'r'i — ^"2) and p = '\J\{pi — P2) and where Hq2 differs from hi + /12 by 

the omission of the center-of-mass HO term of nucleons 1 and 2. Since Vil'"^ acts on 
relative coordinate, the 5'i2 is independent of the two-nucleon center of mass and the 
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two-nucleon center-of-mass Hamiltonian cancels out in Eq. ([27)) . We can see that for 
A > 2 the solutions of ((28|) are bound. The relative-coordinate two-nucleon HO states 
used in the calculation are characterized by quantum numbers \nlsjt) with the radial 
and orbital HO quantum numbers corresponding to coordinate r and momentum p. 
Typically, we solve the two-nucleon Hamiltonian ([28]) for all two-nucleon channels up to 
j = 8. For the channels with higher j only the kinetic-energy term is used in the many- 
nucleon calculation. The model space P2 is defined by the maximal number of allowed 
HO excitations -/Vi2max from the condition 2n + I < -/Vi2max- In order to construct the 
operator oj (PT|) we need to select the set of eigenvectors /C. We select the lowest states 
obtained in each channel. It turns out that these states also have the largest overlap 
with the P2 model space. Their number is given by the number of basis states satisfying 

2n + l < iVi2max- 

An improvement over the two-body effective interaction approximation is the use 
of three-body effective interaction that takes into account the full space three-nucleon 
correlations. If the NNN interaction is included, the three-body effective interaction 
approximation is rather essential for A > 3 systems. First, let us consider the case with 
no NNN interaction. The three-body effective interaction can be calculated as 

T/NN _ 
''3off,123 — 

(29) P3 [e-^^23(/i, + h2 + h3 + V^l'-^ + Vll'^ + V^l^)e^^^''^ - {hi + /i2 + /^a)] ^^3 • 

Here, Sf<^ — arctanh(a-'i23 — '^123) ^^d P3 is a three-nucleon model space projector. 
The P3 space contains all three-nucleon states up to the highest possible three-nucleon 
excitation, which can be found in the P space of the A-nucleon system. For example, for 
A = 6 and A^max — 6 (6?if2) space we have P3 defined by A^i23max — 8. Similarly, for the 
p-shell nuclei with A> 7 and Amax = 6 {QTiH.) space we have Ai23max = 9. The operator 
1^123 is obtained with the help of Eq. (pij) from exact solutions of the Hamiltonian 
hi + h2 + + Vil'"^ + V[l'^ + V^'^, which are found using the antisymmetrized three- 
nucleon Jacobi coordinate HO basis. In practice, we again do not need to calculate ^123, 
rather we apply Eqs. and (pS)) with the three-nucleon solutions. The three-body 
effective interaction is then used in A-nucleon calculations using the effective Hamiltonian 

(30) Hl^^ = + ' 

i=l i<j<k 

where the factor takes care of over-counting the contribution from the two-nucleon 
interaction. 

If the NNN interaction is included, we need to calculate in addition to the 
following effective interaction 



T/NN+NNN _ p 
''3cff,123 ~ ^3 



C.NN + NNN n A O zl (~t A c-NN-f 

{hi +h2 + /^3 + V^12' + ^13 + "^^23' + 14^NN,123)e^-3 



(31) -{hi + h2 + h3)]P3 
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This three-body effective interaction is obtained using fuU space solutions of the Hamil- 
tonian hi + h2 + h3 + V^'^ + V[^'^ + V^'^ + Vnnn,123 • The three-body effective interaction 
contribution from the NNN interaction we then define as 



for,\ T/NNN _ T^NN-l-NNN t/NN 

V'J^J ''3Gff,123 — ''3off,123 ''3cfl,123 • 

The effective Hamiltonian used in the A-nucleon calculation is then 

A A A 

(33) -ffA,off = ^ + ^ _ 2 ^ ^sls.ijk + ^ ^3^S%k ■ 

i—1 i<j<k i<j<k 

At this point we also subtract the i?cM and, if the Slater determinant basis is to be used, 
we add the Lawson projection term P{Hcm — ^hft) to shift the spurious CM excitations. 

It should be noted that all the effective interaction calculations are performed in 
the Jacobi coordinate HO basis. As discussed above, two-body effective interaction is 
performed in the \nlsjt) basis and the three-body effective interaction in the \NiJT) 
basis p3|) . In order to perform the A-nucleon calculation in the Slater determinant HO 
basis as is typically done for A > 4, the effective interaction needs to be transformed to 
single-nucleon HO basis. This is done with help of the HO wave function transformations 
([6]). The details for the three-body case in particular are given in Refs. [23] and [24] , 

As a final remark, we note that the unitary transformation performed on the Hamil- 
tonian should be also applied to other operators that are used to calculate observables. 
If this is done, the model-space-size convergence of observables improves. More details 
on calculation of effective operators are given in Ref. [5S]. 

2'4. Convergence tests. - In this subsection, we give examples of convergence of ab 
initio NCSM calculations. In Fig. [U we show the convergence of the '^H ground-state 
energy with the size of the basis. Thin lines correspond to results obtained with the 
NN interaction only. Thick lines correspond to calculations that also include the NNN 
interaction. The full lines correspond to calculations with two-body effective interaction 
derived from the chiral effective field theory (EFT) NN interaction of Ref. [55] discussed 
in more details in the next section. The dashed lines correspond to calculations with the 
bare, that is the original unrenormalized, chiral EFT NN interaction. The bare NNN 
interaction is added to either the bare NN or to the effective NN interaction in calcu- 
lations depicted by thick lines. Here, we use the chiral EFT NNN interaction that will 
be discussed in details in the next section. We observe that the convergence is faster 
when the two-body effective interaction is used. However, starting at about iVmax — 24 
the convergence is reached also in calculations with the bare NN interaction. The rate 
of convergence also depends on the choice of the HO frequency. In general, it is always 
advantageous to use the effective interaction in order to improve the convergence rate. 
It should be noted that in calculations with the effective interaction, the effective Hamil- 
tonian is different at each point as the effective interaction depends on the size of the 
model space given by iVmax- The calculation with the bare interaction is a variational 
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calculation converging from above with A'max and HO frequency Q as variational param- 
eters. The calculation with the effective interaction is not variational. The convergence 
can be from above, from below or oscillatory. This is because a part of the exact effective 
Hamiltonian is omitted. The calculation without NNN interaction converges to the '^H 
ground-state energy —7.852(5) MeV, well above the experimental —8.482 MeV. Once the 
NNN interaction is added, we obtain —8.473(5) MeV, close to experiment. As discussed 
in the next section, the NNN parameters were tuned to reproduce '^H and '^He binding 
energies. 
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Fig. 1. - ground-state energy dependence on the size of the basis. The HO frequency of 
hQ = 28 MeV was employed. Results with (thick lines) and without (thin lines) the NNN inter- 
action are shown. The full lines correspond to calculations with two-body effective interaction 
derived from the chiral NN interaction, the dashed lines to calculations with the bare chiral NN 
interaction. 



In Fig. [21 we show convergence of the ''He ground-state energy. The NCSM calcula- 
tions are performed in basis spaces up to iVmax = 20. Thin lines correspond to results 
obtained with the NN interaction only, while thick lines correspond to calculations that 
also include the NNN interaction. The dashed lines correspond to results obtained with 
bare interactions. The full lines correspond to results obtained using three-body effective 
interaction (the NCSM three-body cluster approximation). It is apparent that the use 
of the three-body effective interaction improves the convergence rate dramatically. We 
can see that at about Vmax = 18 the bare interaction calculation reaches convergence 
as well. It should be noted, however, that p-shell calculations with the NNN interac- 
tions are presently feasible in model spaces up to Vmax — 6 or iV,nax = 8. The use 
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of the three-body effective interaction is then essential in the p-sheh calculations. The 
calculation without NNN interaction was done for two different HO frequencies, it is 
apparent that convergence to the same result is in both cases. We note that in the case 
of no NNN interaction, we may use just the two-body effective interaction (two-body 
cluster approximation), which is much simpler. The convergence is slower, however, see 
discussion in Ref. [57]. We also note that ''He properties with the chiral EFT NN in- 
teraction that we employ here were calculated using two-body cluster approximation in 
Ref. [28] and present results are in agreement with results found there. Our '*He results 
ground state energy results are —25.39(1) MeV in the NN case and —28.34(2) MeV in 
the NN-I-NNN case. The experimental value is —28.296 MeV. We note that the present 
ab initio NCSM '^H and "'He results obtained with the chiral EFT NN interaction are in 
a perfect agreement with results obtained using the variational calculations in the hyper- 
spherical harmonics basis as well as with the Faddeev-Yakubovsky calculations published 
in Ref. [5S]. A satisfying feature of the present NCSM calculation is the fact that the rate 
of convergence is not affected in any significant way by inclusion of the NNN interaction. 
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Fig. 2. - *He ground-state energy dependence on the size of the basis. The HO frequencies of 
hn = 28 and 36 MeV was employed. Results with (thick lines) and without (thin lines) the 
NNN interaction are shown. The full lines correspond to calculations with three-body effective 
interaction, the dashed lines to calculations with the bare interaction. For further details see 
the text. 



As yet another example of convergence of ab initio NCSM calculations, we present the 
excitation energy calculations of the five lowest excited states of ^Li using the chiral EFT 
NN potential. The NCSM excitation energy dependence on the basis size is presented in 
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Fig. 3. - Calculated positive-parity excitation spectra of ®Li obtained in QhQ,-lAhVl basis spaces 
using two-body effective interactions derived from the chiral EFT NN potential are compared 
to experiment. The HO frequency of TiQ, = 12 MeV was used. 



Fig. [3] for the HO frequency of Ml — \2 MeV. Due to the complexity of the calculations, 
the lowest four state were obtained in basis spaces up to A'max — 14 while we stopped 
at iVmax = 12 for the 2+1 and the 1^0 state. The calculations were performed using the 
two-body effective interaction in the Slater determinant HO basis with the shell model 
code Antoine [30]. Results for other HO frequencies were pubHshed in Ref. ^H- We 
observe that the convergence rate with TVmax is different for different states. In particular, 
the 3+0 state and the 0+1 state converge faster in the higher frequency calculations 
{Ml — 12, 13 MeV), while the higher lying states converge faster in the lower frequency 
calculations {fi^ = 8, 10 MeV). The results in Fig. [3] demonstrate a good convergence of 
the excitation energies in particular for the 3+0 and 0+1 states. An interesting results 
is the overestimation of the 3+0 excitation energy compared to experiment. It turns out 
that this problem is resolved once the NNN interaction is included in the Hamiltonian. 
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3. Light nuclei from chiral EFT interactions 

Interactions among nucleons are governed by quantum chromodynamics (QCD). In 
the low-energy regime relevant to nuclear structure, QCD is non-perturbative, and, there- 
fore, hard to solve. Thus, theory has been forced to resort to models for the interaction, 
which have limited physical basis. New theoretical developments, however, allow us con- 
nect QCD with low-energy nuclear physics. The chiral effective field theory (xEFT) [31] 
provides a promising bridge. Beginning with the pionic or the nucleon-pion system |32| 
one works consistently with systems of increasing nucleon number [33l[34l[35]. One makes 
use of spontaneous breaking of chiral symmetry to systematically expand the strong in- 
teraction in terms of a generic small momentum and takes the explicit breaking of chiral 
symmetry into account by expanding in the pion mass. Thereby, the NN interaction, the 
NNN interaction and also ttN scattering are related to each other. The xEFT predicts, 
along with the NN interaction at the leading order, an NNN interaction at the 3rd order 
(next-to-next-to-leading order or N^LO) [5TI [5^ [57] . and even an NNNN interaction at 
the 4th order (N^LO) [551 ■ The details of QCD dynamics are contained in parameters, 
low-energy constants (LECs), not fixed by the symmetry. These parameters can be con- 
strained by experiment. At present, high-quality NN potentials have been determined 
at order N'^LO [55]. A crucial feature of xEFT is the consistency between the NN, NNN 
and NNNN parts. As a consequence, at N^LO and N^LO, except for two LECs, assigned 
to two NNN diagrams, the potential is fully constrained by the parameters defining the 
NN interaction. 

We adopt the potentials of the xEFT at the orders presently available, the NN at 
N^LO of Ref. [55] and the NNN interaction at N^LO [31137]. Since the NN interaction 
is non-local, the ab initio NCSM is the only approach currently available to solve the 
resulting many-body Schrodinger equation for mid-p-shell nuclei. We are in a position 
to use the ab initio NCSM calculations in two ways. One of them is the determination 
of the LECs assigned to two NNN diagrams that must be determined in A > 3 systems. 
The other is testing predictions of the chiral NN and NNN interactions for light nuclei. 

The NNN interaction at N^LO of the xEFT comprises of three parts: (i) The two- 
pion exchange, (ii) the one-pion exchange plus contact and the three-nucleon contact, 
see Fig. The LECs associated with the two-pion exchange also appear in the NN 
interaction and are therefore determined in the A — 2 system. The one-pion exchange 
plus contact term (D-term) is associated with the LEG cd and the three-nucleon contact 
term (E-term) is associated with the LEC ce- The cd and ce LECs, expected to be of 
order one, can be constrained by the A = 3 binding energy. We then still need additional 
observable to determine the two parameters. Their determination from three-nucleon 
scattering data is difficult due to a correlation of the binding energy and, e.g. the 
nd doublet scattering length [371 and, in general, due to the lack of an in-depth three- 
nucleon scattering phase shift analysis. We therefore investigate sensitivity of the A > 3 
nuclei properties to the variation of the constrained LECs. 

Before presenting results of the ab initio NCSM calculations with the xEFT NN-I-NNN 
interactions, most of which were published in Ref. [39_, let us discuss briefly a few tech- 
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nical details of the calculations with the NNN interaction. The NNN interaction is 
symmetric under permutation of the three micleon indexes. It can be written as a sum 
of three pieces related by particle permutations: 

(34) W = Wi + W2 + W3 

To obtain its matrix element in an antisymmetrized three-nucleon basis we need to con- 
sider just a single term, e.g. Wi. Using the basis introduced in Eq. p^ . a general matrix 
element can be written as 

{NiJT\W\N'i'JT) = ■i{NiJT\Wi\N'i'JT) 

= 3 ^(n/sjt, NCJ\ \NiJT) {n'l's'j't\N'C'J'\\N'i'JT) 

(35) x{{nlsjtMCJ)JT\Wi\{n'l' s' j't\N' C J')JT) . 

We consider just the most trivial part of the xEFT N^LO NNN interaction, the three- 
nucleon contact term and evaluate its matrix element to demonstrate that non-trivial 
effort is needed to include the NNN interactions in many-body calculations. The three- 
nucleon contact term can be written as 

= Et2 ■ f^Sin - r2)S{r3 - n) 
1 1 f 

(36) = Et2 ■ T3j^;^j-^ J d7rid7f2d7r;d7r^|7fi7r2)(7fi7r^| , 

with E = j^j- where is the chiral symmetry breaking scale of the order of the p 
meson mass and F.^ = 92.4 MeV is the weak pion decay constant. The tti and 7r2 are 
Jacobi momenta associated with the the Jacobi coordinates and ^2 defined in Eq. (8) 
and (9). 

The contact term must be regulated before it can be used in many-body calculations. 
We consider a regulator depending on momentum transfer: 

l^™"*.Q = i?f2-T3^^ I d7rid7r2d7rld7ri|7ri7r2)i^(g2;A)F(Q'2;A)(7ri7r^| 
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(37) = Et2 ■ T3 J de;d6|a6>^o(A/2a; A)Zo(|^ei + ^§61; A)(a6| , 

with the regulator function F{q'^;A) = exp(— g*/A"'). We defined momenta transferred 
by nucleon 2 and nucleon 3: Q ^ p'2 ^ P2 = ^ t^i) + '^{'^2 ~ ^2) and Q' — 



— P3 = — Y 3 (■'?2 ~ ''^'i)- Also, we introduced the fimction 
(38) Zo(r; ^) ^ ^ J dqq^joiqr)F{q^; A) . 

This results in an interaction local in coordinate space because of the dependence of the 
regulator function on differences of initial and final Jacobi momenta. We can see that 
after the regulation, the form of the contact interaction becomes much more complicated. 
The three-nucleon matrix element of the regulated term is then obtained in the form 

{{nlsjt,^CJ)JT\W^°'''^\n'l's'j't',Af'C'J')JT) 



E6Sss'ii' +^+ t <j ; 



I ( t t' 1 1 f < t' 1 

T 



1 1 r ^ 1 1 



2 2 2 ) '< 2 2 

xfj'jj'i'c'i-iY- l + J'-J+i + t^ + s 

v I X \ \ J i' X \ \ J' J X 



xY^{-irx'[ 

X ^ 



J ]' s \ \ J' J J \ \ C C \ 

x(roxo|/o)(/:'oxo|£o) 

X j d^id^2ei^lR^i{^i,h)R^c[^2.b)Rn'v{^i,b)RN'C'{^2,b) 
(39) xZo(A/2a;A)^o,x(V^ei, v1'6;A) , 



with a new function 



(40) Zo,x(ri,r2;A) j dqq^Jxiqri)jx{qr2)Fiq^; A) . 



and customary abbreviation I = \^2l + 1. Evaluation of the other N^LO NNN terms is 
still more complicated. 

It is important to note that our NCSM results through A = 4 are fully converged in 
that they are independent of the iVmax cutoff and the hft HO energy. This was demon- 
strated in the subsection on the ab initio NCSM convergence tests in particular for the 
chiral EFT interactions we are investigating here. For heavier systems, we characterize 
the approach to convergence by the dependence of results on A'max and hfl. 

Fig. [5] shows the trajectories of the two LECs cd — ce that are determined from 
fitting the binding energies of the A = 3 systems. Separate curves are shown for 
and ^He fits, as well as their average. There are two points where the binding of ''He 
is reproduced exactly. We observe, however, that in the whole investigated range of 
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Fig. 5. - Relations between cd and ce for which the binding energy of (8.482 MeV) and ^He 
(7.718 MeV) are reproduced, (a) ^He ground-state energy along the averaged curve, (b) ^He 
charge radius Va along the averaged curve. Dotted lines represent the uncertainty due to the 
uncertainties in the proton charge radius. 



CD — Ce, the calculated ^He binding energy is within a few hundred keV of experiment. 
Consequently, the determination of the LECs in this way is likely not very stringent. We 
therefore investigate the sensitivity of the p-shell nuclear properties to the choice of the 
CD — Ce LECs. First, we maintain the A ~ 3 binding energy constraint. Second, we limit 
ourselves to the cd values in the vicinity of the point cu ~ 1 since the values close to the 
point Cd ^ 10 overestimate the *He radius. Also this large value might be considered 
"unnatural" from the xEFT point of view. 

While most of the p-shell nuclear properties, e.g. excitation spectra, are not very 
sensitive to variations of cd in the vicinity of the ^ 1 point, we were able to identify 
several observables that do demonstrate strong dependence on cd- For example, the 
^Li quadrupole moment that changes sign depending on the choice of cd- In Fig. [S] we 
display the ratio of the B(E2) transitions from the ""^"B ground state to the first and the 
second 1+0 state. This ratio changes by several orders of magnitude depending on the 
Cd variation. This is due to the fact that the structure of the two 1+0 states is exchanged 
depending on cd- 

In Fig. [21 we present the ^^C B(M1) transition from the ground state to the 1+1 
state. The B(M1) transition inset illustrates the importance of the NNN interaction in 
reproducing the experimental value f40j. Overall our results show that for < — 2 
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Fig. 6. - Dependence on the cd with the ce constrained by the ^4 = 3 binding energy fit for 
diflferent basis sizes for B(E2;3+0 ^ l+0)/B(E2;3+0 iJO) ratio. The HO frequency of 
HQ = 14 MeV was employed. 



the ^He radius and the ^Li quadrupole moment underestimate experiment while for 
CD > the lowest two 1+ states of are reversed and the ^^C B(Ml;0+0 1+1) is 
overestimated. We therefore select c/j — — 1 as globally the best choice and use it for our 
further investigation. 

We present in Figs. [8] and [9] the excitation spectra of as a function of A^max for 
both the chiral NN-I-NNN, as well as with the chiral NN interaction alone. In both cases, 
the convergence with increasing N^s^x is quite quite reasonable for the low-lying states. 
Similar convergence rates are obtained for our other p— shell nuclei. 

A remarkable feature of the ^"B results is the observation that the chiral NN in- 
teraction alone predicts incorrect ground-state spin of ^"B. The experimental value is 
3+0 while the calculated one is 1+0. On the other hand, once we also include the chiral 
NNN interaction in the Hamiltonian, which is actually required by the xEFT, the correct 
ground-state spin is predicted. Further, once we select the cd value as discussed above, 
i.e. CD = —1, we also obtain the two lowest 1+0 states in the experimental order. 

We display in Fig. [TO] the natural parity excitation spectra of four nuclei in the middle 
of the p— shell with both the NN and the NN+NNN effective interactions from xEFT. 
The results shown are obtained in the largest basis spaces achieved to date for these nuclei 
with the NNN interactions, A^max — 6 (QhH,). Overall, the NNN interaction contributes 
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Fig. 7. - Dependence on the cn with the ce constrained by the ^4 = 3 binding energy fit for 
different basis sizes for the ^^C B(Ml;0+0 1+1). The HO frequency of HQ, = 15 MeV was 
employed. In the inset, the convergence of the B(Ml;0+0 1^1) is presented for calculations 
with (using cd ~ —1) and without the NNN interaction. 



significantly to improve theory in comparison with experiment. This is especially well- 
demonstrated in the odd mass nuclei for the lowest few excited states. The case of 
the ground state spin of ^"B and its sensitivity to the presence of the NNN interaction 
discussed also in Figs. [5] and [S] is clearly evident. We note that the results with the 
NN interaction only in Fig. [5] were obtained with the HO frequency of Sfi = 15 MeV, 
while those in Fig. [10] with hfl = 14 MeV. A weak HO frequency dependence of the 
-^max = 6 results is evident. The results with NN+NNN interaction presented in 
Figs. [9] and [10] were obtained using the same HO frequency. Still, one may notice small 
differences of the Nnmx = 6 results. The reason behind those differences is the use of two 
alternative D-term regularizations (both depending on the momentum transfer, details 
will be discussed elsewhere). As the dependence on the regulator is a higher order effect 
than the xEFT expansion order used to derive the NNN interaction, these differences 
should have only minor effect. It is satisfying that the present -'^"B results appear to 
support this expectation. 

Concerning the ^^C results, there is an initial indication that the chiral NNN interac- 
tion is somewhat over-correcting the inadequacies of the NN interaction since, e.g. 1+0 
and the 4+0 states in ^^C are not only interchanged but they are also spread apart more 
than the experimentally observed separation. In the ^^C results, we can also identify an 



Ab initio NO-CORE SHELL MODEL CALCULATIONS FOR LIGHT NUCLEI 



21 



10 



B 



2^1 
2^ 



3+ 



1+ 

3+ 



NN 
hn =15 MeV 



2^1 
4^ 



3^ 



Exp 



6hQ. 4hQ. 2hQ. 



OhQ. 



Fig. 8. - ^'^B excitation spectra as function of the basis space size A'^max at HQ, — 15 MeV 
using the chiral NN interaction and comparison with experiment. The isospin of the states not 
explicitly depicted is T=0. 



indication of an overly strong correction arising from the chiral NNN interaction as seen 
in the upward shift of the | state. However, the experimental | may have significant 
intruder components and is not well-matched with our state. In addition, convergence for 
some higher lying states is affected by incomplete treatment of clustering in the NCSM. 
This point will be elaborated upon later. These results required substantial computer 
resources. A typical iVmax = 6 spectrum shown in Fig. 1101 and a set of additional exper- 
imental observables, takes 4 hours on 3500 processors of the LLNL's Thunder machine. 
The ^-nucleon calculations were performed in the Slater determinant HO basis using 
the shell model code MFD [JT]. More details on some of the results discussed here are 
published Ref. [39]. 

The calculations presented in this section demonstrate that the chiral NNN inter- 
action makes substantial contributions to improving the spectra and other observables. 
However, there is room for further improvement in comparison with experiment. In these 
calculations we used a strength of the 27r-exchange piece of the NNN interaction, which is 
consistent with the NN interaction that we employed (i.e. from Ref. 26 ). This strength 
is somewhat uncertain (see e.g. Ref. [H]). Therefore, it will be important to study the 
sensitivity of our results with respect to this strength. Further on, it will be interesting 
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Fig. 9. - ^'^B excitation spectra as function of the basis space size TVmax at hO, — 14 MeV using 
the chiral NN+NNN interaction and comparison with experiment. The isospin of the states not 
explicitly depicted is T=0. 



to incorporate sub-leading NNN interactions and also four-nucleon interactions, which 
are also order N-^LO Finally, it will be useful to extend the basis spaces to iVmax = 8 
(8hfl) for A > 6 to further improve convergence. 

4. Cluster overlap functions and S-factors of capture reactions 

In the ab initio NCSM calculations, we are able to obtain wave functions of low-lying 
states of light nuclei in large model spaces. An interesting and important question is, 
what is the cluster structure of these wave functions. That is we want to understand, how 
much, e.g. an ^Li eigenstate looks like ''He plus deuteron, an ^Be eigenstate looks like 
^He plus '^He, an ®B eigenstate looks like ''Be plus proton and so on. This information 
is important for the description of low-energy nuclear reactions. To gain insight, one 
introduces channel cluster form factors (or overlap integrals, overlap functions). The 
formalism for calculating the channel cluster form factors from the NCSM wave functions 
was developed in Ref. |42j. Here we just briefly repeat a part of the formalism relevant 
to the simplest case when the lighter of the two clusters is a single-nucleon. 

We consider a composite system of A nucleons, i.e. *B, a nucleon projectile, here a 
proton, and an 1-nucleon target, i.e. ''Be. Both nuclei are assumed to be described by 
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Fig. 10. - States dominated by p-shell configurations for ^"B, ^^B, ^"^C, and ^"'C calculated at 
iVmax = 6 using hO. = 15 MeV (14 MeV for ^°B). Most of the eigenstates are isospin T=0 or 
1/2, the isospin label is explicitly shown only for states with T=l or 3/2. The excitation energy 
scales are in MeV. 



eigenstates of the NCSM effective Hamiltonians expanded in the HO basis "with identical 
HO frequency and the same (for the eigenstates of the same parity) or differing by 
one unit of the HO excitation (for the eigenstates of opposite parity) definitions of the 
model space. The target and the composite system is assumed to be described by "wave 
functions expanded in Slater determinant single-particle HO basis (that is obtained from 
a calculation using a shell model code like Antoine). 
Let us introduce a projectile-target "wave function 

(6 . . .ei-2r'f|<i>j;\-;;_;;);'*^;<5,) = Y.[jmI,M,\JM){lmilm,\jm)'-^ 

(41) y-Yi^,[r)XmA^i ■ ■■U-2\A - lahAh) , 

"where (^i . . . ^a-2|^ — la/iMi) and Xm^ ^-re the target and the nucleon "wave function, 
respectively. Here, I is the channel relative orbital angular momentum, ^ are the target 
Jacobi coordinates defined in Eq. ([7]) and r — (ri + r2 + • • • + M-i) — de- 

scribes the relative distance bet"ween the nucleon and the center of mass of the target. 
The spin and isospin coordinates "were omitted for simplicity. 
The channel cluster form factor is then defined by 

(42) 9tt^,.A-^.i.(^ = iA^J\^\!^:X'^^' 

"with A the antisymmetrizer and |ylAJ) an eigenstate of the A- nucleon composite sys- 
tem (here ^B). It can be calculated from the NCSM eigenstates obtained in the Slater- 
determinant basis from a reduced matrix element of the creation operator. The derivation 
is as follo-ws. First, we use the relation p?|) for both the composite A-nucleon and the 
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target A — 1-nucleon eigenstate. With the help of relations analogous to ([6]): 
y^^{LMlm\Qq)(pNLMiRcM) {ta) = 

Mm 

(43) {n'l'N'L'Q\NLnlQ)_^il'm'L'M'\Qq)ipr,,i,m'iU-i)VN'L'M'i^o) , 

n'l'm' N' L' M' 

we obtain 

(44) SD {AX J\A^\f,^'^:lY^ ; nl)su = {nmi\00nll) _^ {A\J\A^\f,^]:l\-l ^ ' 

with a general HO bracket due to the CM motion. The nl in ([44]) refers to a replacement 
of Sr by the HO i?„i(r) radial wave function. Second, we relate the SD overlap to a 
linear combination of matrix elements of a creation operator between the target and the 
composite eigenstates sBiAXJlal^^jlA — la/i)sD- The subscript SD refers to the fact 
that these states were obtained in the Slater determinant basis. Such matrix elements 
are easily calculated by shell model codes. The result is 

^ ' {i^,j);ali ^ '''^ ' {nl00l\00nll)_L_ ' 

(45) X sD{A^J\\aiij\\A-laIi)sD . 

The eigenstates expanded in the Slater determinant basis contain CM components. A 
general HO bracket, which value is simply given by 

(46) {nmi\00nll) _^ = (-1)' 

then appears in Eq. in order to remove these components. The Rni{r) in Eq. 

is the radial HO wave function with the oscillator length parameter b = ^ where 

m is the nucleon mass. 

A conventional spectroscopic factor is obtained by integrating the square of the cluster 
form factor; 

(47) = / drA4lU-i.i^^'')\" ■ 

A generalization for projectiles {— the lighter of the two clusters) with 2, 3 or 4 
nucleons is straightforward, although the expressions become more involved. In all cases, 
the projectile is described by wave function expanded in Jacobi coordinate HO basis, 
while the composite and the target eigenstates are expanded in the Slater determinant 
HO basis. Full details are given in Ref. [42] . 

As an example, in Fig. [11] we present the cluster overlap function of ""^^B ground state 
with ^Li+^He. Results are given for ^Li in the 1+0 ground state and in the 3+0 excited 
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state. We can see that the ^Li 1+0 ground-state component is rather small. The 
ground state is dominated by a superposition of S-, D- and G- waves of relative motion 
of *He and ^Li in the 3"'"0 state. 

The overlap functions introduced in this subsection are relevant for description of 
low-energy 7-capture reactions important for nuclear astrophysics. Next, we investigate 
three reactions of this type. 

4'1. Be(p,j)^ B . - The ^Be(p,7)®B capture reaction serves as an important input 
for understanding the solar neutrino flux ^43]. Recent experiments have determined the 
neutrino flux emitted from with a precision of 9% [33]. On the other hand, theoretical 
predictions have uncertainties of the order of 20% [33 135] . The theoretical neutrino flux 
depends on the ''Be(p,7)^B S-factor. Many experimental and theoretical investigations 
studied this reaction. 

In this subsection, we discuss a calculation of the ^Be(p,7)®B S-factor starting from 
ab initio wave functions of ^B and ^Be. It should be noted that the aim of ab initio 
approaches is to predict correctly absolute cross sections (S-factors), not only relative 
cross sections. The full details of our ^Be(p,7)*B investigation were published in Refs. [47l 
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EH]. 

Our calculations for both ''Be and nuclei were performed using the high-precision 
CD-Bonn 2000 NN potential [49] in model spaces up to lOhil (A^max = 10) for a wide 
range of HO frequencies. From the obtained ®B and ^Be wave functions, we calculate 
the channel cluster form factors (overlap functions, overlap integrals) g^i'i'lj.^_i^j_^{r) 
as discussed in the previous subsection. Here, A = 8, Z is the channel relative orbital 
angular momentum and r = (ri + ^2 + • • ■ + ^^4-1) ~ t^a describes the relative 

distance between the proton and the center of mass of ^Be. The two most important 
channels are the p- waves, 1 = 1, with the proton in the j = 3/2 and j = 1/2 states, 
j = I + s, s = 1/2. In these channels, we obtain the spectroscopic factors of 0.96 and 
0.10, respectively. The dominant j = 3/2 overlap integral is presented in Fig. [T2l by 
the full line. The 10?iJ7 model space and the HO frequency of hft = 12 MeV were used. 
Despite the fact, that a very large basis was employed in the present calculation, it is 
apparent that the overlap function is nearly zero at about 10 fm. This is a consequence 
of the HO basis asymptotic behavior. As already discussed, in the ab initio NCSM, 
the short-range correlations are taken into account by means of the effective interaction. 
The medium-range correlations are then included by using a large, multi-Tiil HO basis. 
The long-range behavior is not treated correctly, however. The proton capture on ^Be 
to the weakly bound ground state of ®B associated dominantly by the El radiation is a 
peripheral process. In order to calculate the S-factor of this process we need to go beyond 
the ab initio NCSM as done up to this point. We expect, however, that the interior part 
of the overlap function is realistic. It is then straightforward to find a quick fix and 
correct the asymptotic behavior of the overlap functions, which should be proportional 
to the Whittaker function. 

One possibility we explored utilizes solutions of a Woods-Saxon (WS) potential. In 
particular, we performed a least-square fit of a WS potential solution to the interior 
of the NCSM overlap in the range of — 4 fm. The WS potential parameters were 
varied in the fit under the constraint that the experimental separation energy of ^Be-|-p, 
Eq = 0.137 MeV, was reproduced. In this way we obtain a perfect fit to the interior of 
the overlap integral and a correct asymptotic behavior at the same time. The result is 
shown in Fig. [13 by the dashed line. 

Another possibility is a direct matching of logarithmic derivatives of the NCSM over- 
lap integral and the Whittaker function: ■^ln{rgij{r)) = ■jpln{CijW^ri.i+i/2{'2kor)) , 
where 77 is the Sommerfeld parameter, ko = \/2^Eo/h with /i the reduced mass and 
Eq the separation energy. Since asymptotic normalization constant (ANC) Cij cancels 
out, there is a unique solution at r = Rm- For the discussed overlap presented in Fig.lT^ 
we found Rm = 4.05 fm. The corrected overlap using the Whittaker function matching is 
shown in Fig. [T^] by a dotted line. In general, we observe that the approach using the WS 
fit leads to deviations from the original NCSM overlap starting at a smaller radius. In 
addition, the WS solution fit introduces an intermediate range from about 4 fm to about 
6 fm, where the corrected overlap deviates from both the original NCSM overlap and the 
Whittaker function. Perhaps, this is a more realistic approach compared to the direct 
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Whittaker function matching. In any case, by considering the two alternative procedures 
we are in a better position to estimate uncertainties in our S-factor results. 
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Fig. 12. - Overlap function, rg{r), for the ground state of with the ground state of '^Be plus 
proton as a dependence on separation between the ^Be and the proton. The p-wave channel 
with j = 3/2 is shown. The full line represents the NCSM result obtained using the CD- 
Bonn 2000 NN potential, the lOhQ model space and the HO frequency of hQ = 12 MeV. 
The dashed lines represent corrected overlaps obtained from a Woods-Saxon potential whose 
parameters were fit to the NCSM overlaps up to 4.0 fm under the constraint to reproduce the 
experimental separation energy. The dotted lines represent overlap corrections by the direct 
Whittaker function matching. 



In the end, we re-scale the corrected overlap functions to preserve the original NCSM 
spectroscopic factors (Table 2 of Ref. [47|). In general, we observe a faster convergence 
of the spectroscopic factors than that of the overlap functions. The corrected overlap 
function should represent the infinite space result. By re-scaling a corrected overlap 
function obtained at a finite iVmax, we approach faster the infinite space result. At the 
same time, by re-scaling we preserve the spectroscopic factor sum rules. 

The S-factor for the reaction ^Be(p,7)*B also depends on the continuum wave func- 
tion, i?y ^ . As we have not yet developed an extension of the NCSM to describe continuum 

wave functions (see, however, the discussion in Sect. [5]), we obtain i?!^' for s and d waves 
from a WS potential model. Since the largest part of the integrand stays outside the 
nuclear interior, one expects that the continuum wave functions are well described in 
this way. In order to have the same scattering wave function in all the calculations, we 
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chose a WS potential from Ref. [50j that was fitted to reproduce the p-wave l"*" reso- 
nance in ®B. It was argued [51] that such a potential is also suitable for the description 
of s- and d-waves. We note that the S-factor is very weakly dependent on the choice of 
the scattering-state potential (using our fitted potential for the scattering state instead 
changes the S-factor by less than 1.5 eV b at 1.6 MeV with no change at MeV). 

Our obtained S-factor is presented in Figs. [T3] where contribution from the two partial 
waves are shown together with the total result. It is interesting to note a good agreement 
of our calculated S-factor with the recent Seattle direct measurement [35]. 
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Fig. 13. - The ■'Be(p,7)'*B S-factor obtained using the NCSM overlap functions with corrected 
asymptotics as described in the text. The dashed and dashed-dotted lines show the contribution 
due to the i = 1, j = 3/2 and j = 1/2 partial waves, respectively. Experimental values are from 
Refs. miMl. 



In order to judge the convergence of our S-factor calculation, we performed a detailed 
investigation of the model-space-size and the HO frequency dependencies. We used the 
HO frequencies in the range from Ml = 11 MeV to Ml ~ 15 MeV and the model 
spaces from QMl to \QMl. By analysing these results, we arrived at the S-factor value of 
S'i7(10keV) = 22.1 ± 1.0 eV b. 

4'2. ^ He(a,^ f Be. - The '^He(Q;,7)^Be capture reaction cross section was identified 
the most important uncertainty in the solar model predictions of the neutrino fluxes in 
the p-p chain [46]. We investigated the bound states of ^Be, ■^He and ^He within the ah 
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initio NCSM and calculated the overlap functions of ^Be bound states with the ground 
states of ^He plus ''He as a function of separation between the '^He and the a particle. 
The obtained p-wave overlap functions of the ^Be 3/2" ground state excited state are 
presented in Fig.[T3]by the full line. The dashed lines show the corrected overlap function 
obtained by the least-square fits of the WS parameters done in the same way as in the 
^B^^Be+p case. The corresponding NCSM spectroscopic factors obtained using the 
CD-Bonn 2000 in the lOhV. model space for ^Be {Uhn for ^•'^Re) and HO frequency of 
hfl = 13 MeV are 0.93 and 0.91 for the ground state and the first excited state of ^Be, 
respectively. We note that contrary to the *B^''Be-|-p case, the '^Be^^He-|-a p-wave 
overlap functions have a node. 
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Fig. 14. - The overlap function, rgir), for the first excited state of ^Be with the ground state 
of ^He plus a as a dependence on separation between the ^He and the a particle. The p-wave 
channel overlap function with j — 3/2 is shown. The full line represents the NCSM result 
obtained using the CD-Bonn 2000 NN potential and the lOhfl model space for ^Be {I2hfl for 
^'*He) with the HO frequency of hQ = 13 MeV. The dashed line represents a corrected overlap 
obtained with a Woods-Saxon potential whose parameters were fit to the NCSM overlap up to 
3.4 fm under the constraint to reproduce the experimental separation energy. 

Using the corrected overlap functions and a "^He-I-Q! scattering state obtained using 
the potential model of Ref. [54 we calculated the •^He(Q;,7)^Be S-factor. Our lOhil 
result is presented in the left panel of Fig. [151 We show the total S-factor as well as 
the contributions from the capture to the ground state and the first excited state of 
^Be. By investigating the model space dependence for Shfl and lOTifl spaces we estimate 
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the ■^He(a,7)^Be S-factor at zero energy to be higher than 0.44 keV b, the value that we 
obtained in the discussed case shown in Fig. [151 Our results are similar to those obtained 
by K. Nollett [55] using the variational Monte Carlo wave functions for the bound states 
and potential model wave functions for the scattering state. 
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Fig. 15. - The full line shows the ^He(a,7)^Be S-factor obtained using the NCSM overlap 
functions with corrected asymptotics. The dashed lines show the ^Be ground- and the first 
excited state contributions. The calculation was done using the CD-Bonn 2000 NN potential 
and the lOhCl model space for ^Be {UhO. for ^'''He) with the HO frequency of HQ = 13 MeV. 



4'3. ^H(a,^)''Li. - An important check on the consistency of the '^He(a,7)^Be S- 
factor calculation is the investigation of the mirror reaction ■^H(a,7)^Li, for which more 
accurate data exist [56]. Our results obtained using the CD-Bonn 2000 NN potential 
are shown in Fig. [Tni It is apparent that our ^H(a,7)^Li results are consistent with our 
^He(a,7)''Be calculation. We are on the lower side of the data and we find an increase 
of the S-factor as we increase the size of our basis. 

More details on the ah initio NCSM investigation of the ■^He(a,7)^Be and ■^H(Q;,7)''Li 
S-factors are given in Ref. [57] . 

5. Towards the ah initio NCSM with continuum 

In the previous section, we highlighted shortcomings of the ah initio NCSM, its incor- 
rect description of long-range correlations and its lack of coupling to continuum. If we 
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Fig. 16. - The full line shows the ^H(Q:,7)^Li S-factor obtained using the NCSM overlap functions 
with corrected asymptotics. The dashed lines show the '^Li ground- and the first excited state 
contributions. The calculation was done using the CD-Bonn 2000 NN potential and the lORfi 
model space for ^Li {12hn for and *He) with the HO frequency of = 13 MeV. 



want to build upon the ab initio NCSM to microscopically describe loosely bound sys- 
tems as well as nuclear reactions, the approach must be augmented by explicitly including 
cluster states such as, e.g. those given in Eq. (|¥T1) . and solve for their relative motion 
while imposing the proper boundary conditions. This can be done by extending the ab 
initio NCSM HO basis through the addition of the cluster states. This would result in an 
over-complete basis with the cluster relative motion wave functions as amplitudes that 
need to be determined. The first step in this direction is to consider the cluster basis 
alone. This approach is very much in the spirit of the resonating group method (RGM) 
[S5] , a technique that considers clusters with fixed internal degrees of freedom, treats the 
Pauli principle exactly and solves the many-body problem by determining the relative 
motion between the various clusters. In our approach, we use the ab initio NCSM wave 
functions for the clusters involved and the ab initio NCSM effective interactions derived 
from realistic NN (and eventually also from NNN) potentials. 

The general outline of the formalism is as follows. The many-body wave function is 
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approximated by a superposition of binary cluster channel wave functions 
with 

(49) $1^'^'") = V^^'^Vt^<5(r - TA-a^a) ■ 

Here, A is the antisymmetrizer accounting for the exchanges of nucleons between the two 
clusters (which are already antisymmetric with respect to exchanges of internal nucleons). 
The relative-motion wave functions ipi, depend on the relative-distance between the center 
of masses of the two clusters in channel v. They can be determined by solving the many- 
body Schrodinger equation in the Hilbert space spanned by the basis functions (|49p : 



(50) h¥^^ = 



(A) 



E 



df 



■7^(^-«-«)(r r ) - EU^J-^'^Hr', r)] Mr) 



where the Hamiltonian and norm kernels are defined as 

(51) nli-''^''Hf',r) = ($;;l7"'"^ \aha 

(52) M^-^^^Hr',r) = U[iV-^ A' 



^(A-a,a) 



The most challenging task is to evaluate the Hamiltonian kernel and the norm kernel. 
We now briefly outline, how this is done when ab initio NCSM wave functions are used 
for the binary cluster states. From now on, let us consider the cluster states with a 
single- nucleon projectile (a = 1 in Eq. |48)) . A generalization is straightforward. Using 
an alternative coupling scheme compared to Eq. (|4ip . we introduce 



(A-1,1)JMTMt , 



= ^ (/i Ml i m, I sm) (smZm; I JM) (Ti Mt, i mt I TMt) ' ^ 
(53) X iiA-i)Xm^Xm, (Ci . . . U-2\A - lahMiTiMr,) , 

with the spin and isospin coordinates omitted to simplify the notation. The Jacobi 
coordinates were defined in Eq. ([7]). Using the latter cluster basis and the following 
definition of the antisymmetrizer A — 1/VA{1 — ^f^i Pj,A) with Pj^a the transposition 
operator of nucleons j and A, the norm kernel can be expressed as 



S{r' 



{A-l)J2Rn'i'{r') 



(54) 



X 11. nn^l I \PA,A^l\<f) 



{aliTu 



2 2 



)sl 



;nl)Rni{r) 
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with /i = (a'/(T{, ii)s', = {aliTi, ^^)s and Pa^a-i the transposition operator of 



nucleons A and A — 1. The coordinates r are related to £,a~i by r y ^:j£,A-i 
and the HO length parameter of the radial HO wave functions is 6 = ^ ^-'^m^ - '^^^'^ 
matrix element of the transposition operator Pa,a~-i can be directly evaluated using 
the ab initio NCSM wave functions expanded in Jacobi coordinate HO basis following 
a procedure analogous to the derivation of Eq. (fH|) . However, a crucial feature of the 
ab initio NCSM approach is that the matrix elements that enter the norm kernel and 
the Hamiltonian kernel can be equivalently evaluated using the ab initio NCSM wave 
functions expanded in the Slater determinant HO basis. This is achieved in two stages. 
First, we calculate the SD matrix element as 

X ss'jj'iff(-i)^;+-'"+''(-i)^i+ 5 +^ 

(55) X SD (A - la'/; T{ 1 1 1 (a^ 5^ i ) (^^) \\\A- lahT^)sB ■ 

Second, it is possible to show that the matrix element in the SD basis is related to the 
one in the Jacobi coordinate basis: 



(56) 



^ _ 11, ,-,n'l'\PA.A^M^''''^i\, ,;nl)sB^ 



2 2 



\" -1-1, 22'" 2 2'" 

E/^iA-l,l)J^T ,,, IP ,^(A-1,1)J,^T 

nArTi'l'Jr -^1^1^ 2 2'''5 V (q/iTi, 2 jjS'r 

^ ;7' 72/ 1^s+^,.-s-^' J '/r 1 J s' Jr \ 

X 2^ii j,(-ij 1 L J z M L J r r 

NL ^ ^ ^ ' 

X {nXN Ll\QQnll) _!^{n%N Ll'\QQn'l'l') . 



This relation then defines a matrix that one inverts to get the Jacobi-coordinate matrix 
element. This is analogous to what was done to obtain the translationally invariant 
density in Ref. |59| . The Hamiltonian kernel can be evaluated in a similar yet more 
involved way. It consists of a kinetic term, a NN potential direct term associated with 
the operator Va,a-i[^ — Pa,a-i) and a NN potential exchange term associated with 
the operator Va,a-2Pa,a-i (plus terms arising from the NNN interaction). The ability 
to employ wave functions expanded in the SD basis opens the possibility to apply this 
formalism for nuclei with A > 5. 

In Fig. I17i we show the exchange part of the norm kernel for the n+*He system, in 
particular the second term of Eq. (|54p multiplied by rr' . It is apparent that we are able 
to reach convergence for the kernel. Furthemore, the ^<S'i/2 channel shows, as it should 
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Fig. 17. - The exchange part of the norm kernel of the n+^'He system. Left, the convergence 
with the size of the basis of the *He wave function for the ■^S'1/2 channel. Right, results for 
channels are compared. The chiral EFT NN potential was used. 



the effect of the Pauli principle. Indeed, the '^He w&Ye function is dominated by the four 
nucleon s-shell configuration. The Paufi principle prevents adding the fifth nucleon to 
the same shell. 

In Fig. [181 we show the direct and the exchange contributions of the NN potential 
to the Hamiltonian kernel as well as their sum for the n+^He system. Again, the Pauli 
principle is manifest in the ^6*1/2 channel. We were able to obtain the presented results 
using wave functions expanded both in the Jacobi-coordinate and the SD basis. The two 
independent calculations gave identical results as expected. A converged calculation of 
the ^51/2 phase shift together with experimental data is presented in Fig. [191 Full details 
regarding this approach are given in Ref. [5Tj . 

6. Conclusions 

The ah initio NCSM evolved into a powerful many-body technique. Presently, it is 
the only method capable to use interactions derived within the chiral EFT for systems of 
more than four nucleons, in particular for mid-p-shell nuclei. Among its successes is the 
demonstration of importance of the NNN interaction for nuclear structure. Applications 
to nuclear reactions with a proper treatment of long-range properties are under develop- 
ment. Extension to heavier nuclei is achieved through the importance-truncated NCSM 
|62j . Within this approach, ab initio calculations for nuclei as heavy as ""^Ca become 
possible. 

* * * 

I would like to thank all the collaborators that contributed to the cited papers and. 
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Fig. 18. - The NN potential direct and exchange terms of the Hamiltonian kernel of the n-|-*He 
system. The and ^P3/2 channels are compared. The low- momentum NN potential of 

Ref. [BD] was used in the 14hQ model space with the HO frequency of HQ. = 18 MeV. 
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